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Abstract 

We present a calculation of the effective weak Hamiltonian which 
governs 6 — > s 7 and b ^ s g transitions in two different renormal- 
ization schemes (NDR and HV). In the leading logarithmic approxi- 
mation, we show that the coefficients of the effective Hamiltonian are 
scheme independent only when one takes correctly into account the 
scheme dependence of one- and two- loop diagrams. We demonstrate 
that in NDR there are contributions which were missed in previous 
calculations. These contributions are necessary to obtain scheme in- 
dependent coefficients in the final results. 
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1 Introduction 



Radiative decays of B mesons have been the subject of several investigations 
in the last years With the exception of two recent papers by M. 

Misiakp, |^ and |T0[ however, none of these studies has taken into account the 
full set of leading logarithmic corrections. Moreover calculations made in the 
dimensional reduction scheme (DRED)[0, ^ Q seem to differ from the results 
obtained in the naive dimensional regularization scheme (NDR)[|l], ^, ^, ^]. 

In refs.0 and the calculation of refs.[|I], ^] was repeated in NDR and 
found in agreement with the original calculation. However the origin of the 
difference between the results obtained in NDR and DRED has never been 
clarified. In the literature it is always implicitely (or explicitely) assumed 
that the anomalous dimension matrix which enters in b radiative decays is 
regularization scheme independent, as usually happens in leading logarithmic 
calculations. Thus, bearing in mind possible errors in the calculation, the 
difference between NDR and DRED has remained so far a mistery. 

In this paper we have computed the anomalous dimension matrix relevant 
for 6 — i> s 7 and b ^ s g decays in two different regularization schemes: 
NDR and the t'Hooft-Veltman scheme (HV)[0. We show that the values 
of some diagrams are indeed regularization scheme dependent. We explain 
the origin of the scheme dependence and show that the regularization scheme 
independent effective Hamiltonian is obtained only by taking into account the 
scheme dependence of some one loop diagrams, which appear at zero order in 
as- The contributions of these diagrams were included in the calculation of 
Misiakp, 1^, but their role was not fully understood. We show that the easiest 
way to obtain the scheme independent result is to use the HV regularization 
scheme. We also briefly discuss the role of "effervescent" operators and of 
operators which vanish by the equations of motion. Finally we add some 
comments on the gauge invariance of the final result. 
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2 Evolution and Initial Conditions for the Co- 
efficients 



The effective Hamiltonian for 6 — s 7, (7 decays can be written as: 



where Vij are the elements of the CKM|T^, |T3[ mixing matrix. The operator 
basis is given by: 

Ql = {SaCi3){V-A){Ci3ba){V-A) 
Q2 = iSaCa){V-A)iCi3bi3)(^V-A) 

^3,5 = {Saba){V-A) i^pQfs) (VtA) 

q=u,d,s,--- 

a){V^A) 

q=u,d,s,--- 

^ Qd(^ - U TP 

~ YQ^'i^b^ (^''^ {V +A) 

Qs = Y^^bSaCT'^v+A)ta/3b^G^u (2) 

In eq.(0) the subscript (V^A) indicates the chiral structure, a, [3 are colour 
indices, denotes the b quark mass, = — | is the electric charge of the 
down-type quarks and g (e) is the strong (electro-magnetic) coupling. Our 
colour matrices are normalized in such a way that Tr{t^t^) = 6^^ /2. 

The coefficients C{fi) obey the renormalization group equations: 

where t = ln{My^/ and as = g'^/'^n. The factor of 2 in eq.(^ normal- 
izes the anomalous dimension matrix as in refs.|I^-|T^. 7-^ includes the 
contribution due to the renormalization of mj, and, in the case of Qs, the 
contribution due to the renormalization of the gluon field and of the strong 
couphng constant g, see e.g. ref. |]T8|, ^. 

The anomalous dimension matrix is "almost" triangular, in the sense 
that the operators of dimension six can mix with the magnetic operators Qr 
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and Qs, whereas the dimension five operators Qy g can only mix between 
themselves. For later convenience we introduce a reduced 6x6 anomalous 
dimension matrix, which mixes the operators Qi,Q2, ■ ■ ■ jQe among them- 
selves and a 6-component column vector Cr{fJ.) = (Ci(//), C2(/x), . . . ,Cq{ij,)) 
on which 7^ acts. We also introduce two 6-component vectors, related 
to the two-loop anomalous dimension matrix, (3^ = (717, 727, • • • , 76?) and 

P8 = (718,728, • • • ,768)- 

In terms of these quantities the renormalization group equations can be 
written as: 



d 



(4) 

where = q;s(/x) and n'^d/diJ? = ii'^d/dn^ + (3{as)d/das. One can easily di- 
agonalize the sub-matrix of magnetic operators and write the renormalization 
group equations as: 



d 



as 
An 
a. 



2n'^-^V8{i^) = ^788^;8(a*) 



dfj,^ 



47r 



(5) 



where: 



with 



0:7 



Criii) + dr ■ Criii) + 

— * 

Csifj) + ds • Crin) 
07 + 



787 



777 - 788 



(6) 



= (777! - 7r) 

= (7881 - 7^) ^ 0s 



787 



777 - 788 



(7) 
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Figure 1: One loop diagrams which mix Q5 and Qe with Qj and Qs- The mix- 
ing appears at order and the coefficients of the mixing are regularization 
scheme dependent. 

Contrary to the claims of refs. and Ps do depend on the 

regularization, as our calculations in NDR and HV explicitly show, see be- 
low. A regularization dependent anomalous dimension matrix in a leading 
logarithmic (but two loop) calculation is an exotic phenomenon and deserves 
some explanations. At one loop, but at order the operators Q5 and Qq 
can mix with Qj and Qs with finite coefficients, through the diagrams shown 
in fig.|l]. In the figure, the cross denotes a mass {rrib) insertion in the heavy 
quark propagator of the penguin loop. The coefficients of the mixing are reg- 
ularization dependent. They vanish in HV (and DRED), but not in NDR. As 
a consequence, the values of all two loop diagrams contributing to Py and Ps, 
which have as sub-diagrams those in fig.|l], will depend on the regularization. 
Furthermore the diagrams in fig.^ which vanish in HV, have instead to be 
considered in NDR. 

An example of a diagram which depends on the regularization is reported 
in fig.|^. The integration over the penguin loop, where a bare four fermion 
operator is inserted, will transform it in an effective S'-y^D'^t^G^^.b local vertex 
(Irl = I'^i^ ^ 75))- This effective vertex is, by the equations of motion, 

equivalent to s'y'^t^b J2q \ (jlltiLt^Q + QI/iRt^o) , which can be written as a 
linear combination of Q3, . . . ,Qq. Thus, the divergent part of the coefficient 
of the effective vertex combines with the finite, regularization dependent, 
mixing coefficient of Qs.e with Qr^s- This means that it will give a scheme 
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Figure 2: Two loop diagram contributing to the mixing of the operators 
Qifi with and Qs- The value of the diagram is regularization scheme 
dependent. This figure explains the scheme dependence of the diagram. 
The integration over the penguin loop mixes the original operator Qi with 
s^'^D^G'^j^t^b, denoted by a big square. By the equations of motion, this 
operator is equivalent to a four fermion operator, denoted by a small square 
in the figure, see also the text. 

dependent contribution to and Ps- 

We now show how, from the finite one loop coefficients and the two 
loop anomalous dimension matrix, which are separately scheme dependent, 
a physical, scheme independent, effective Hamiltonian is obtained. Let us 
consider the finite one loop mixing of the operators Q5 and Qq with Q^: 

< s^\Heff\b > = Cy{h) < s-f\Q7{n)\b > +C5{n) < s^\Q5{n)\b > 

+CM<sj\QM\b> (8) 

We can interpret the non vanishing of the one loop matrix elements of 
and Qq as the effect of a mixing matrix acting at order q;°, see below: 

< s^\Heff\b >= drill) < s^\Q7{ii)\b > (9) 
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where Cji^fi) = Cj{n) + Z-j ■ Cr^jj). The vector Zj (and Zg, see below) is 
regularization dependent. Z-j vanishes in HV (and DRED) and 



Z7 = (0,0,0,0,2,2iV) 



(10) 



in NDR. Eq. (^) corresponds to a finite renormahzation of the operators such 
that the matrix elements of and Qq vanish. In this way the renormalized 
operators are the same in NDR and HV. Similarly one finds: 

<sg\Heff\h> = Csifi) < sg\Qs{fi)\b > +C5ifJ.) < sg\Q5{fi)\b > 

= CM < sg\QM\h> (11) 

with CM = C8(/i) + Zs ■ Crifi). ^8 = in HV (and DRED) and 

fg^ (0,0,0,0,2,0) (12) 

in NDR. 

The renormahzation group equations (^ can be easily solved and com- 
bined with eqs. (|^) and ( p!TD to obtain the effective Hamiltonian expressed in 
terms of operators renormalized at the scale /i: 



H^ff ~ g^(/i)-a.(^) 

787 



7r - 777 1 



/?7 + 



787 



+ 

788 - 777 



111 - 788 
^8(yU) + (ir - 788 i) /38 " C7r(/i) 



/?8 ) + ^7 



d(/i) 



7r- 788l) /38 + ^8 



787 



788 - 777 



-^8(/i) + 



7r - 777 1 



-1 



Pi 



787 



777 - 788 



(13) 



+ 



787 



777 - 788 



a(/i) 



787 



788 - 777 
'^8 + ^8 



7r - 788 Ij Ps + Zs 

d(/i)|Q8(/x) 



QM 
(14) 



Cr{^), v-j{n) and v%{n) are scheme independent quantities because they 
obey the scheme independent renormahzation group equations (|^). /Jy^s and 
^7^8 are scheme dependent quantities. We now prove that the combinations 



787 



^ ^ V 777 - 788 / 777 - 788 

^^8 = (ir - 7881) + ^8 (15) 



appearing in eq. ([T^) are instead scheme independent. Let us introduce the 
renormahzation matrix Z"" which gives the renormalized operators Q{^) in 
terms of the bare operators Qs, in the regularization labeUed as "a" : 

Q{li) = (z^y'QB (16) 



The anomalous dimension matrix is defined in terms of Z°- as: 

T = 2{z\^^)y\'-^^z\^,) (17) 

If one uses two different regularization schemes, "a" and "6" say, the 
renormahzation constants are related through the equation: 

Z"" = Z^f (18) 

from eqs.(0) and (p!8| ) we obtain: 

In other for example for the operators appearing in the AS* = 1 

weak Hamiltonian, f differs from the identity matrix by terms of order as, 
i.e. f = 1 + fracas4:nAr |T^-|T^. In the present calculation instead, because 
of the diagrams in fig.|l|, f differs from 1 by terms of order see eqs.(^ and 



TTl). The matrix f can be written as: 

_ A 7 \ 

(20) 



h -AZ 

12 



where 16,2 are 6x6 and 2x2 identity matrices and AZ is a 6 x 2 matrix. 
Its two columns are given by the difference of the vectors Z-i and Zg in the 
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regularizations "a" and "6", AZ = (Z^ - Z^, Z^ - Z\). At the order we are 
interested in, eq. (|T9D gives: 



7 -7 



(21) 



Eg. (|2T|) , which connects the anomalous dimension matrices in two different 
regularizations, holds diagram by diagram and can be used as a check of 
the calculation of any single diagram in different regularization schemes [|17 



More details will be given in ref.[21 



Taking into account the structure of the matrix AZ we find: 



t = % 



(A/37), = A7,7 
(A/^s), = A7,8 



777I - 7rj AZ7 + 787AZ8 
7881 - 7^) Af8 



(22) 



where j denotes the components of the vectors A/37,85 J = 1? • • • ? 6. Eq.(|2^) 
demonstrates that the combinations Q-j and appearing in eq. (|T5|) are regu- 
larization scheme independent. Since Z-j and Z^ are zero in HV, in the present 

and 



NDR 
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HV 



case AZ is simply given by Z7 8 computed in NDR (A7 = 7 
AZ = Z^^^ — Z^^ = Z^^^). We have computed the anomalous dimension 
in both NDR and HV 0. From the one- loop diagrams in fig.|I| and the two- loop 
diagrams in fig.^ and ^, we obtain the anomalous dimension matrices (P3|), 
(|26D, given below. Combining (|25D and ([26| ) with eqs.(|l^,|l^), one indeed 
finds the same u;7,8 iii the NDR and HV schemes. 

In refs.[0, |], ^, they considered a restricted operator basis, in the NDR 
scheme. Our results for the two loop Feynman diagrams agree with those 
reported there. The anomalous dimension matrix given in these references is 
however regularization scheme dependent. The scheme dependence follows 
from the non- vanishing of the diagrams in fig.|l|, which spoil the equations of 
motion of the regularized theory: 



(23) 



■^The reduced mixing matrix 7^ has been recently computed up to the next-to-leading 
order in refs-lUl-lfiTl. 
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Figure 3: Diagrams in the reduced (a) and complete (6) bases. In (a) we 
insert the operator s'-f^D'^t'^G^j^b, in (6) the corresponding four fermion op- 
erators. Because of mass insertions, which spoil the equations of motion, in 
NDR the two diagrams have different values. On the contrary, they give the 
same result in the HV scheme. 

Indeed, using the NDR scheme, the values of the two diagrams in fig.^ differ, 
because of the contributions coming from the insertion of in the loop 
propagators, see also fig.0. It is easy to show that the scheme dependence 
can be eliminated by enforcing the equations of motion with a suitable, finite 
subtraction, at order 

In NDR, when one uses the complete basis (0), the diagrams in fig.^ 
have to be computed[§]. These diagrams vanish in the HV scheme. In NDR 
instead they give a non-zero contribution to the anomalous dimension, pre- 
cisely for the same reason that the diagrams in fig.|l] do not vanish. In the 
final answer ([T^), the contribution of each of the diagrams in fig.^ is exactly 
cancelled by the corresponding term in Zj and Zs, as it can be checked dia- 
gram by diagram. In fact the diagrams in fig.^give a contribution which goes 
like — (^7^ — 777,881) Zj^g. The only exceptions are the diagrams P4 and F4 

in fig.^ (see also fig.0). In this case, Zj^g cancel the regularization dependent 
part of that two-loop diagram, giving the same result one would obtain in 

^This diagram would be present also in the reduced operator basis used in refs.||l|, |[ ||]. 
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Figure 4: Diagrams which contribute both in the NDR and HV schemes. In 
the evaluation of the diagrams counter-terms and "effervescent" operators 
have been included, see also ref. pll . 



HV. In HV, since Z-j and .Zg vanish, one obtains the regularization scheme in- 
dependent result directly from the evaluation of the diagrams in fig.§. Those 
calculations, || which did not include the effect of § in the effective 
Hamiltonian, are not correct because of the regularization dependence of the 
diagrams P4 and F4 in fig.^. On the other hand, for the complete basis, it is 
necessary to include the contribution of the diagrams in figj^. 

We now summarize our discussion. In HV, where Z7 § vanish, we have 
only to compute the contribution to the anomalous dimension matrix of the 
diagrams in fig.§ (besides the usual one loop diagrams which mix the opera- 
tors Qi, . . . , Qe)- The resulting anomalous dimension matrix is regularization 
independent and does not require any further manipulation. In NDR, the 
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two loop anomalous dimension has to be computed from the diagrams in 
figs. I and ^ The result is then combined with Zj^s as in eq.(14). The final 
result is then identical to the result obtained in HV, as we have explicitely 
checked. This implies that our results for the anomalous dimension matrices 
in NDR and HV, (p5D-(^), satisfy the general relation (|2T|) . We report below 
the anomalous dimension matrix in HV, to be used for the evolution of the 



coefficients, together with the initial conditions [|20| 



Ci, C3, 



,C,{Mw) 
Cs{Mw) 



where x = mf/M^. 



1 




_ 2a;2 

'2(1 



Inx 



2fl 



X 



Inx + 



Sx'^ + 5x — 7x 
4(1 

x^ — 5x'^ — 2x 
4(1 - x)3 



(24) 
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+ 



2nf_ 
3 



12 



la 



47r 



8 

9 2N 



12Qi, jV^-i 
<3d 2iV 



22jV 
9 



58 
9Af 



232 
9 2Ar 



8^^J^ Af2_i 12n^(Af^-l) 
9 2Ar 2Ar 



_ 116 _|_ fi„ 



12 + 



9 m 

22Nn 



L 



9N 



-16 



2W 



8n/ iv2_i _ 12nj(JV^-l) 
~9~ 2Ar 2iV 

oiV2-l 



V 



' 2N 

' 2N 



-47V + A _ 6n, 



32Nnf , 50nf 
9 9iV 





AN- 



(25) 



where = + is the number of active flavours, A*" is the number of 
colours, a = 1, . . . , 8, i = 1, . . . , 6 and a — 7,S. In the NDR scheme we 
obtain 7^°^ = 7"^ and: 

I ai I at 

/ 6 



7 



as 
An 





-16 _|_ 12Qu N^-l 

9 2N ^ Qd 2N 

184 
9 2N 



l&Uf Af2-l I 12n/(jV^-l) 
'~9 2N ' 2N 



22Ar 
9 



46 
9Ar 



12 + 



_ 92. , c 
9 9iV ^ 

22Nnf _ m-ri; 



40 



Ar2-i 



9 



9N 



16n/ N^—l 
~9 2W 



2iV 

12n/(Af^-l) _|_ 40Af(Af^-l) 



V 



2N 
' 2iV 
' 2Af 



2W 



4Ar - ^ - 6ny 



o _ 32Nnf , 62n/ 
"° 9 ~'~ 9iV 



4Ar 



where nj = rid 



Qd 
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Figure 5: Diagrams which have to be considered only in the NDR scheme. 
In the evaluation of the diagrams counter-terms and "effervescent" operators 
have been included. 

Let us compare our results with previous calculations. 

In refs. it was assumed that the values of the two loop diagrams 

were scheme independent. In the NDR scheme, we agree with the results of 
refs.[|], H, 1^ when we reduce to the basis used in those references. However 
the necessity of adding Zy g, to obtain a scheme independent result, was 
missed. 

In ref.]^ Misiak realized the necessity of adding Z^, even if he did not 
understand its role for the scheme independence of the final result. Moreover 
the paper contains several errors: counter-terms of current-current diagrams 
and "effervescent" operators^ were not taken into account and there are in- 

'^Without effervescent operators refation (|2l| ) between HV and NDR would fail. 
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consistencies between different anomalous matrix elements, which can be 
related to the same set of diagrams. In ref. 0, some of these errors have 
been corrected: counter-terms of current-current diagrams and effervescent 
operators have been considered and included in the final result. However, 
our results (p6[) and the matrix given in eq.(21) of ref. do not agree. By 
looking to colour factors, the disagreement seems to originate from diagrams 
P2 and P3 in fig.H with a mass {mi,) insertion in the loop propagators. 

The calculation done in DRED in ref. 0, ^ ^ could have lead in principle 
to the correct results. The authors however were biased by the idea that 
the two loop calculation had to be regularization scheme independent and 
interpreted the difference between NDR and DRED as a failure of the DRED 
scheme. We have demonstrated that there are not priviliged schemes and 
that all the schemes, if correctly used, give the same results for the Wilson 
coefficient^. We are actually performing the calculation also in the DRED 
scheme to show that this is indeed the case. 

We were unable to understand how the calculation of ref. [0 was done. 



We disagree with the final result for the anomalous dimension matrix given 



in this reference. The results of ref. 10 differ also from the results of ref. 



We want to add a few comments on the operator basis and gauge in- 
variance. In refs. and a redundant basis of operators which vanish (or 
coincide) by the equations of motion was used. In ref. [|r^ , it was shown that 



this leads to an unnecessary complication and that the correct result can be 
obtained by working directly with a reduced basis of indipendent operators, 
as the basis in eq.(^. This remains true in the present calculation, as we 
have explicitly checked. As noticed already in ref. P, P] , it is not necessary to 
work in the background field gauge to obtain a gauge invariant result. We 
have done (for the non-abelian diagrams) the calculation in the Feynman and 
background gauges (see also ref. (jT^) with identical final results. We plan to 
present a phenomenological analysis in a separate publication, together with 
all the details of the calculation (treatment of the effervescent operators, 
counter-terms, gauge invariance, contribution of any single diagram, etc.) 
and the results in the DRED scheme. 



"At this order in the perturbative expansion. 
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